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We develop expressions for the nonlinear wave damping and frequency correction of a field of 
random, spatially homogeneous, acoustic waves. The implications for the nature of the equilibrium 
spectral energy distribution are discussed. 



On 

On 



3 



> ■ 

o : 
o ■ 
r» ■ 
o : 

On ■ 

: 

>v 
■ 

o : 

o ■ 



INTRODUCTION AND GENERAL 
DISCUSSION 



Weak or wave turbulence, which describes the behav- 
ior of a spatially homogeneous field of weakly interact- 
ing, random dispersive waves, has led to spectacular suc- 
cess in our understanding of spectral energy transfer pro- 
cesses in plasmas, oceans and planetary atmospheres [0. 
Furthermore, the subject provides a useful paradigm for 
helping one think about some of the challenges of fully de- 
veloped turbulence. First and foremost, the equation for 
the long time behavior of the spectral cumulants (which 
are in one to one correspondence with the spectral mo- 
ments) are closed without making a priori and unjustifi- 
able assumptions on the statistics of the processes (such 
as the quasi-Gaussian approximation). Second, the ki- 
netic equation, which describes the spectral energy trans- 
fer via n-wave resonant processes, admits classes of exact 
equilibrium solutions that can be identified as pure Kol- 
mogorov spectra, namely equilibria for which there is a 
constant spectral flux of one of the conserved densities 
(e.g. energy, number density) of the physical process 
under consideration. By contrast, the thermodynamic 
equilibria, which have very little relevance in any turbu- 
lence theory that must account for a sink at small scales, 
have zero flux. Third, the theory allows one to glimpse 
the origin of the intermitency and the breakdown of the 
conditions under which one can expect relaxation to one 
of the finite flux Kolmogorov equilibria. 

The basic ideas for writing down the kinetic equation 
to describe how weakly interacting waves share their en- 
ergies go back to Peicrls but the modern theories have 
their origin in the works of Hasselman Q, Benney and 
Saffmann fl, Kadomtsev Zakharov [y, Benney and 
Newell A particularly important event in this his- 

tory was the discovery of the pure Kolmogorov solution 
by Zakharov Usually, the thermodynamic equilib- 
rium solutions can be seen from the kinetic equation by 
inspection. On the other hand, the solutions, correspond- 
ing to pure Kolmogorov spectra are much more subtle 
and only emerge after one has exploited scaling symme- 
tries of the dispersion relation and the coupling coefficient 
via what is now called the Zakharov transformation JjJ . 

But success to this point, namely the natural closure, 



depended crucially on the fact the waves were dispersive. 
This means that the group velocity is neither constant 
in amplitude nor direction, or alternatively stated, the 
dispersion tensor 



Dap = 



dk a dkp 



< a, (3 < a 



(1.1) 



has full rank. Here d is the system dimension, Greek let- 
ters (here a and (3) denote tensor indices varying from 1 
to space dimension d, and 



: W (k) 



(1.2) 



is the linear dispersion relation. The reason for closure 
is slaving. In a field of weakly interacting random dis- 
persive waves, the first step to slaving is achieved by the 
linear characteristics of the wave trains. Systems, which 
initially are strongly correlated, are decorrelated because 
different waves carry statistically independent informa- 
tion from long distances at different speeds and direc- 
tions. Cumulants of order N > 2 tend to zero on the fast 
time scale ui^ 1 (ui is a typical frequency at which the en- 
ergy is injected). The system approaches a state of exact 
joint Gaussian statistics. The energy at each wavevec- 
tor remains constant and there is no transfer. But the 
systems of interest to us are nonlinear and therefore, al- 
though the cumulants undergo an initial decay, they are 
regenerated by the nonlinear terms. In particular, the 
cumulant of the order N is regenerated both by cumu- 
lants of higher order and by and products of lower order 
cumulants. The second important reason for closure is 
the following. The important terms in the regeneration 
of the A^ th order cumulant are not the terms containing 
cumulants of order higher then TV but rather those terms 
which are products of lower order cumulants. Important 
means that even though the nonlinear coupling is weak, 
the effects of these terms persist over long times because 
of resonant interaction. Furthermore, the regeneration 
process takes place on a much longer time scale than does 
the initial decorrelation process due to wave dispersion. 
On this long time scale, namely the time taken for triad 
or quartet (or, as in some rare cases, quintic) resonances 
to produce order one effects, the system of equations for 
the cumulant hierarchy becomes closed. If e is a typical 
dimensionless wave amplitude (for acoustic waves it is 
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5p/po, the ratio of average fluctuation density amplitude 
to the ambient value) then this time (measured in units 
of the timescale lo^ 1 ) is e~ 2 for triad resonances and e~ 4 
for quartet resonances, although there is an additional 
frequency correction in the latter case that comes in on 



the 



time scale. 



Mathematically, these results are obtained by pertur- 
bation theory, in which the terms leading to long time cu- 
mulative effects can be identified, tabulated and summed. 
The method closely parallels that of the Dyson- Wyld di- 
agrammatic approach which will be discussed in Section 
4. A key part of the analysis is the asymptotic (limt_>oo) 
evaluation of certain integrals such as 



/ 



/(k r )A 



N 



s r o>(k r ) 



JY 



5 M Udk r 



where 



Jo 



ih 



(1.3) 



(1.4) 



and 5(x) is the Dirac delta function. The function A(h) 
contains the fast (oscillations of the order of ujJ 1 ) time 
i, whereas the other functions in the integrand, here de- 
noted by /(k r ), only change over much longer times. The 
exponent of A(h) is X^i^ Srw(k r ) where ui(k r ) is the lin- 
ear dispersion relation and s r (often s r = ±1) denotes its 
multiplicity. For example, in acoustic waves, a wavevec- 
tor k has two frequencies corresponding to waves running 
parallel and anti-parallel to k . The maximum contribu- 
tion to integrals such as (F3) in the limit of large time t 
occurs on the so called resonant manifold M, where 



N 



N 



k r = , h = s r w(k r .) 



(1.5) 



for some choices of the sequence s r . However, the precise 
form of the asymptotic limit also depends on whether the 
zeros of h on M are simple or of higher order. For the 
case of (fully) dispersive waves, such as gravity waves on 
deep water, Rossby waves, waves of diffraction on optical 
beams, the zero of h is simple and (for sufficiently smooth 
/) one has 



hm r fjhf^-'dh 
= 7rsgn(i)/(0) + iP / ^dh 



(1.6) 



or schematically, 



A(h) 0C7rsgn(tW/i) + iP(-), (1.7) 

h 

where P denotes Cauchy principal value. In these cases, 
the integrand in the kinetic equation, the equation de- 
scribing the resonant transfer of spectral density, con- 



tains products of energy densities and the Dirac delta- 
functions <5(X)i Sj-w(k r )) and S(^2^ k r ) clearly indicat- 
ing that spectral energy transfer takes place on the res- 
onant manifold M. The asymptotic equations for the 
change of the higher order cumulants can be interpreted 
as a complex frequency modification whose real part de- 
scribes the expected nonlinear shift in frequency and 
whose imaginary part describes a broadening of the res- 
onant manifold along its normal directions. 

But acoustic waves are not fully dispersive. The linear 
dispersion relation, 



w(k) = c\k\ = c*lk 2 



k = (fc||, fcj 



(1.8) 



where c is the sound speed, leads to a dispersion tensor 
which has rank (d — 1). As we will see, this changes the 
asymptotic. Furthermore, three wave resonances occur 
between wavevectors which are purely collinear. There- 
fore, since the kinetic equation (KE) only considers wave 
interaction on the resonant manifold, there is no way of 
redistributing energy out of a given direction. At best, 
the KE will only describe spectral energy transfer along 
rays in wavevector space. Moreover, depending on di- 
mension d, the long time behavior of the integrals (1.4) 
differ greatly. For a given vector k, the locus of the reso- 
nant partners k! and k — k x in a resonant triad is given 
by the surface in ki space defined by 



/i(ki) = sifci + s 2 \k - fci| - s\k\ = 



(1.9) 



Here s,Si,S2 = ±1. For d = 1 and the appropriate 
choices of the wave directions si, S2 and s, this manifold 
is all ki. Therefore the fast oscillations in the integral 
are of no consequence and do not cause any decorrelation 
to occur. All waves moving in the same directions travel 
with the same speed. Initial correlations are completely 
preserved. Moreover, we know that for one dimensional 
compressible flow, nonlinear terms, no matter how weak 
initially, eventually lead to finite time multivalued solu- 
tions. Assuming the usual viscous regularization, multi- 
valued solutions are replaced by shocks, namely almost 
discontinuous solutions where discontinuities are resolved 
across very thin viscous layers. One would naturally ex- 
pect an energy spectrum E\{k) which reflects this fact, 
namely 



E 1 (k) tx 1/k 2 



(1.10) 



In two dimensions, one has dispersion (diffraction) in one 
direction. Indeed, for d > 1, while 



V kl h = 



(1.11) 



on the manifold M, the Hessian of h(ki) i s n ot identi- 
cally zero. In two dimensions, the integral ( |l.4| ) behaves 
as 



/(x) 6X ^ ^ X ^ ^-dxcc2t I f(x) e^p(ix 2 t)dx 1 



(1.12) 
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which grows as t 1 / 2 as t — > oo. In three dimensions, the 
growth is much weaker. Since that is the case we will 
look at in detail, we give the exact result. Let 

k=(K> 0,0,0), kj = {K x ,K y ,K z ), 
k 2 = {K-K x ,-K y ,-K z ) . 

Then, for si = s 2 = s, 



h = c(s 1 \k 1 \ + s 2 \k-k 1 \-sK) 



(1.13) 
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(K 2 + K 2 ) + 0(K*K 2 K z ,...) 



2K X {K-K X 
near the resonant value (K, 0, 0). The integral 
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tends to 



af f(K x ,0,0;- S) s)(-K x )(K-K x )dK x 

+ a [ f(K x , 0, 0; s, s)K x (K - K x )dK x 
Jo 

+ a f{K x ,0,0;s,-s)K x (K x -K)dK x (1.14) 

JK 

2ias f° 

log* / f(K x , 0,0,; -s,s)K x (K-K x )dK x 

2ias f K 

log*/ f(K x ,0,0;s,s)K x (K-K x )dK x 
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in the limit t — > oo. Here a = ir 2 /Kc and we have kept 
only the leading order real and im aginar y co ntributions. 
The essential difference with ([li]) and (1.7) is the ad- 
ditional Dirac delta function multiplied by log* in the 
imaginary term. This will not change the kinetic equa- 
tion for the spectral energy density. If we write the total 
energy per unit volume E as 



E = 2p Q c 2 e 2 J e(k)dk 



(1.15) 



where po is the ambient density and e a measure of am- 
plitude, then 



rfe(k) 
dt 



St(e, e) 



St(e,e)=^ +1 ^ K \-2 
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- ( 7 + l)e(k)e( 7 k) 



dka(l — a) e(afc)e((l — a)k) 



;(k)e((l - a)k) - (1 - a)e(fc)e(ak) 



(1.16) 



where p is the adiabatic constant [p — Poip/po)^] an d 
|k| = K . In d dimensions a little calculation show, that 

3-d 



the RHS of ( 1.16 ) has the t dependence t 2 so that in 
general the nonlinear interaction time t^l for the res- 
onant exchange of spectral energy is £ 2 t~ = 0(1) or 

4 

tnl oc (T^=d _ (Note that for d > 5, there is no cumula- 
tive effect of this resonance.) 



While the extra term in (1.14) proportional to zlogt 
plays no role in the spectral energy transfer, it will, how- 
ever, appear in the frequency modification. Calculat- 
ing the long time behavior of the higher order cumulants 
leads to a natural re-normalization of the frequency, 



;(k) = c|k| 



1 - 2tt(p+ l) 2 e 2 ln 



1 



(3 2 e((3k)dl3 



+ 0(e 2 ) 



+ in 2 (p + l) 2 e 2 



P 2 e(pk)dp 



+ TTT /' l f3 3 e{Pk)d[3+ \k\ /' (3e(f3k)df3 



(1.17) 



where k = k/K. The calculation of the frequency re- 
normalization is the new result of this paper. We present 
two derivations of this result, in the framework of the 
above analysis and making use of a diagrammatic per- 
turbation approach. 



The equation ( 1.16 ) is nothing but a "regular" kinetic 
equation for the three-wave interactions, written in a dis- 
persionless limit u — c|k|. In this case three wave reso- 
nant conditions 



±w(k) = ±w(ki)±w(k 2 ) k = ki+k 2 



(1.18) 



can be satisfied if and only if all three vectors k, kx , k 2 are 
parallel, as a result, the integration over ki,k 2 is along 
line parallel to k. It is unclear a priori that the three wave 
kinetic equation can be used in the dispersionless case; is 
certainly less plausible in the two dimensional case where 
the formal implementation of the kinetic equation leads 
to stronger divergences. 

The derivation presented above is taken from the arti- 
cle of Newell and Aucoin j9) , who made the first serious 
attempt of an analytical description of the dispersionless 
acoustic turbulence. 

Newell and Aucoin ||] also argued that a natural 
asymptotic closure also obtains in two dimensions be- 
cause of the relative higher asymptotic growth rates of 
terms in the kinetic equation involving only the spectral 
energy, but this is still a point of dispute, is not yet re- 
solved and will not be addressed further here. 



Independently the kinetic equation (1.16) was applied 
to acoustic turbulence by Zakharov and Sagdeev || who 
used it just as a plausible hypothesis. However, Zakharov 
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and Sagdeev also suggested an explicit expression for the 
spectrum of acoustic turbulence 



e(fc) cx k~ 3/2 



(1.19) 



which is just a Kolmogorov-type spectrum, first obtained 
by Kolmogorov from dimensional considerations in the 
conte xt of hydrodynamic turbulence. Here, however, the 
(1.19) is an exact solution of the equation 



St(e,e) = 0. 



(1.20) 



The proof of this fact can be found in the [gj . One should 
also mention, that the quantum kinetic equation applied 
to a description of a system of weakly interacting dis- 
persionless phonons were done as long ago as in 1937 by 
Landau and Rumer ( |Tc||). 

Kadomtsev and Petviashvili jll[] criticized this result 
on the grounds that the kinetic equation is in the disper- 
sionless case can hardly be justified because of the special 
nature of the linear dispersion relation. They suggested 
that acoustic turbulence in two and three dimensions was 
much more likely to have parallels with its analogue in 
one dimension. We have already mentioned in that case 
that the usual statistical description is inadequate both 
because there is no decorrelation dynamics and because 
shocks form no matter how weak the nonlinearity ini- 
tially is. The equilibrium statistics relevant in that case 
is much more likely to be a random distribution of dis- 
continuities in the density and velocity fields which lead 
to an energy distribution of (1.1C). Further, Kadomt- 
sev and Petviashvili argued that even in two and three 
dimensions one would expect the same result, namely 



fc d-1 e(k) oc k- 



(1.21) 



a random distribution of statistically independent shocks 
propagating in all directions. 

But wave packets traveling in almost parallel direc- 
tions are not independent. Consider a solid angle con- 
taining N = (k\\/k±) d ~ 1 wavepackets with wavevectors 
(fc||,fcjj where fey = l^ 1 is a typical length scale of 
the fluctuating field in the direction of the propagation, 
and fc_L <C fen. The shock time r s h for a single wave 
packet would be l^pN/E cx {I / ce)N i - 1 / 2 \ where E is 
the total energy in the field. The dispersion (diffraction) 
time Tdisp, namely the time over which several different 
packets have time to interact linearly, is of the order of 
fcii/(cfcjj cx IN 2 /^ 1 } /c. As we have already observed, 
the nonlinear resonance interaction time tnl for spec- 
tral energy transfer is {l/c)e~^^ b ~ d ^ . The ration is 
r disp : r sh : r NL = N 2 ^-V : ATV2 £ -i . e -4/( 6 -«0. In 
the limits N — > oo, e — > 0, the shock time is sandwiched 
between the linear dispersion time and nonlinear inter- 
action time and, if we choose N(e) by equating the first 
two, all three are the same. Moreover, the phase mixing 
which occurs due to the crossing of acoustic wave beams, 
occurs on a shorter time scale, a fact that suggests that 
the resonant exchange of energy is the more important 



process. But even then, several very important questions 
remain. 

1 . To what distribution does the energy along a given 
wavevector ray relax? 

2. How does energy become shared between neighbor- 
ing rays? 

3. Does energy tend to diffuse away from the ray with 
maximum energy or can it focus onto that ray? In 
the latter case, one might argue that shock forma- 
tion may again become the relevant process espe- 
cially if the energy should condense on rays with 
very different directions. 

The aim of this paper is to take a very modest first step 
in the direction of answering these questions. In partic- 
ular, we present a curious result. The fact that there is 
a strong (e 2 In 1/e 2 ) correction to the frequency leads us 
to ask if that terms could provide the dispersion required 
to allow the usual triad resonance process carry energy 
between neighboring rays. At first sight, it would appear 
that that is indeed the case, that the modified nonlinear 
dispersion law is 



«(k)=c(k)(l 



e 2 ln\n(k)) 



(1.22) 



where fl is proportional to |k|. But a surprising and 
nontrivial cancellation occurs which means that the first 
corrections to the wave speed still keeps the system non 
dispersive in the propagation direction. 

While this fact is the principal new result of this pa- 
per, our approach lays the foundation for a systematic 
evaluation of the contribution to energy exchange that 
occurs at higher order. Indeed, we expect that some of 
the terms found by Benney and Newell || involving gra- 
dients across resonant manifolds which, in the fully dis- 
persive case, are not relevant because the resonant three 
wave interaction gives rise to an isotropic distribution, 
may be more important in this context. 

The paper is written as follows. In the next Sec- 
tion, we derive the equation of motion for acoustic waves 
of small but finite amplitude. A second approach dis- 
cussed in Subsect. 2B starts from the Hamiltonian for- 
mulation of the Euler equations and again makes use of 
the small amplitude parameter of the problem to simplify 
the interaction Hamiltonian. As we will see in Subsect. 
2C both approaches are equivalent and which approach 
to use is the question of taste. 

Next, in Section 3 we write down the hierarchy of equa- 
tions for the spectral cumulants and solve them pertur- 
batively. Certain resonances manifest themselves as alge- 
braic and logarithmic time growth in the formal pertur- 
bation expansions and mean that these expansions are 
not uniformly asymptotic in time. The kinetic equation, 
describing the long time behavior of the zeroth order 
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spectral energy, and the equations describing the long 
time behavior of the zeroth order higher cumulants are 
simply conditions that effectively sum the effect of the 
unbounded growth terms. Under this renormalization, 
the perturbation series becomes asymptotically uniform. 
By asymptotically uniform, we mean that the asymptotic 
expansion for each of the cumulants remains an asymp- 
totic expansion over long times. All unbounded growths 
are removed. While this procedure in principle requires 
one to identify and calculate unbounded terms to all or- 
ders, in practice one gains a very good approximation by 
demanding uniform asymptotic behavior only to that or- 
der in the coupling coefficient where the unboundedness 
first appears. 

In other words this means that if one finds that if 
the first two terms of the asymptotic expansion are 
1 + e 2 tipi +..., then the effective removal of ipi will remove 
all terms which are powers of (e 2 t) in the full expansion. 
Likewise, it also assumes that there appear no worse sec- 
ular terms at a higher order, such as for example e 4 t 3 ip2- 
To achieve uniformity, one requires an intimate knowl- 
edge of how unbounded growth appears. This sort of per- 
turbative analysis was first done in the thirties by Dyson. 
A technical innovation was to use graph notations, called 
diagrams, for representing lengthy analytical expressions 
for high order terms in the perturbation series. It of- 
ten happens that one can find the principal subsequence 
of terms just by looking on the topological structure of 
corresponding diagrams. This method of treating pertur- 
bation approaches is called the diagrammatic technique. 

The first variant of diagrammatic technique for non- 
equilibrium processes was suggested by Wyld in 
the context of the Naiver Stokes equation for an incom- 
pressible fluid. This technique was later generalized by 
Martin, Siggia and Rose who demonstrated that 

it may be used to investigate the fluctuation effects in 
the low-frequency dynamics of any condensed matter sys- 
tem. In fact this technique is also a classical limit of the 
Keldysh diagrammatic technique H] which is applicable 
to any physical system described by interacting Fermi 
and Bose fields. Zakharov and L'vov |l5) extended the 
Wyld technique to the statistical description of Hamilto- 
nian nonlinear- wave fields, including hydrodynamic tur- 
bulence in the Clebsch variables p6| . In section 4, we 
will use this particular method for treating acoustic tur- 
bulence. 

Note that in such a formulation, unbounded growths 
appear as divergences (or almost divergences) due to 
the presence of zero denominators caused by resonances, 
the very same resonances, in fact, that give rise to un- 
bounded growth in our more straightforward perturba- 
tion approach. Moreover, diagrammatic techniques are 
schematic methods for identifying all problem terms and 
for adding them up. If one uses the diagram technique 
only to the first order at which the first divergences ap- 
pear, this is called the one- loop approximation and is 
equivalent to identifying the first long time nonlinear ef- 
fects. This is exactly analogous to what we will do in 



our first approach in this paper although we will also 
display the diagram technique. The one loop approxima- 
tion will give the same long time behavior of the system 
for times of t^l defined earlier. In Appendix C we an- 
alyze two loop diagrams and show that some of them 
gives the same order contribution to 7^ as two loop dia- 
grams. Nevertheless one may believe, that even one-loop 
approximation gives qualitatively correct description of 
the dynamics of the system. 

The last Section 5 is devoted to some concluding re- 
marks and the identification of the remaining challenges. 
We now begin with deriving the basic equations of mo- 
tion for weak acoustic turbulence. 



2. BASIC EQUATION OF MOTION FOR WEAK 
ACOUSTIC TURBULENCE 



A. Straightforward Approach 



Consider the Euler equations for a compressible fluid: 



dp/dt + V • (pv) = , 

dw/dt + (v • V)v = -Vp{p)/p . 



(2.1) 



Here t>(x, t) is the Euler fluid velocity, p(x,t) the den- 
sity, and p(r, t) is the pressure which, in the general 
case, is a function of fluid density and specific entropy 
s [p = p(p, s)]. In ideal fluids where there is no viscosity 
and heat exchange, the entropy per unit volume is carried 
by the fluid, i.e. obeys the equation ds/dt+ (v • V)s = 0. 
A fluid in which the specific entropy is constant through- 
out the volume is called barotropic; the pressure in such 
fluid is a single- valued function of density p = p(p). 
In this case, Vp/p may be expressed via the gradi- 
ent of specific enthalpy of unit mass w = E + pV and 
dw = Vdp = dp/ p. Thus, Vp/p = Vro. 

Writing the fluid density p(x, t) as po(l + r/(x, t)), the 
velocity field as u(x, t), the pressure field asp = po{l+i])^ 
and the enthalpy as 



dp 



"■ I - = -^ T (l + 0*-l)»7+ 
p fi-1 



(p-l)(p-2) 
o V +•■•) 



one can write ( |2.l| ) to third order in amplitude in the 
following form 



di] dvi d 



dt dxi 
dv 



3X., 



_J_ +C 2^V = _ y dl >* c (M~ 2 ) 9 2 

dt dx m dx m 2 dxi 



c 2 (p-2)(p-3) d ?3 

6 8X n 



(2.2) 



(2.3) 



Let us introduce new variables as 
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r)(x,t)= /^ea s (k,i)e lkx+l ^ (k)t dk, (2.4) 

•* s 

«,(x,i) = JY, ^ea^Kty^^dk, (2.5) 



where < e -C 1, a; (A;) = c|k [ and V connotes summa- 
tion over s = ±1. From ( [2.2| ) and (2.3), 



da s (k,t) 
Ft 



x 6(kp + k g - k) exp{j[s p o;(kp) 4- s q u)(k q ) - sw(k)]t} 



x a s<! (k«j,t)a' !! ''(k r ,t)(5(kp + k g + k r - k) 

x exp{i[s p aj(kp) + s q to(k q ) + s r uj(k r ) — suj(k)]t} (2.6) 

where the summation is done over all signs of s p , s q , s r 
and we used the shorthand notation uj p — Lu(k p ). The 
coupling coefficients are, 



zc f kk i 1 kki-i 



T",c p ,a q / ~P , 

K,K p ,K, 4 \ SpUJp SqLOq SpUJpSqUlq 



- k>, k. 



k,k P ,k„k. = ^(M-2)(M-3) 



(2.7) 
(2.8) 



These coefficients have the following important proper- 
ties: 



• £ k 'k'k is symmetric under the interchange of p 
and q 



9 J^ S P' S ' S q 

kp ,k, kg 



(SpLUp/ SU))L 



k.kp ,kg 



On the resonant manifold M, given by 



i 



-h = s p \kp\ + s q \k - kp| - sjkj = , (2.9) 



(2.10) 



B. Hamiltonian Description of Acoustic Turbulence 

1. Equations of Motion and Canonical Variables 

Co nside r again the Euler equations for a compressible 
fluid (2.1). The enthalpy of a unit mass w = E + pV is 
equal to the derivative of internal energy of unit volume 
e(p) = Ep with respect to fluid density: w — 5e/Sp . As 
a result of direct differentiation with respect to time, it is 
readily evident that equations ( |2.lD conserve the energy 
of the fluid 

ft = J [pv 2 /2 + £(p)]dr . (2.11) 



One can show (and see for example ||1[|) that Eqs. (2.1) 
may be written in the Hamiltonian form: 

dp/dt = 5H/6$ , d$/dt = -5H/5p , (2.12) 
dX/dt = SH/Sp , dp/dt = -5H/8X , (2.13) 



if the velocity v(r, t) is presented in terms of two pairs of 
Clebsch variables (p, $) and (A, v) as follows, 



v = A 



Vp 



(2.14) 



Here the ener gy (|2.1l|) is expressed in terms (p, $) and 
(A, v) so that (2.14 ) bec omes the Hamiltonian of the sys- 
tem. As seen from (2.14), the case with A = or p =const 
corresponds to potential fluid motions which are d efine d 
by a pair of variables (p, $) according to equations (2.12). 
It is convenient to transform in the k-representation from 
the real canonical variables, $(k),/?(k) to the complex 
ones b(k) and 6*(k), 



$(k) = - iygy j[i(k) - 6*(-k)] , 
<5p(k) = v / (M72)[6(k) + 6*(-k)] . 



(2.15) 
(2.16) 



Here Sp(k) — [p(k) — po(k)] is the Fourier transform of 
density deviation from the steady state. 



2. Hamiltonian of Acoustic Turbulence 



where |k| = K. Note that if k = (if, 0,0), the res- 
onant manifold is not of codimension one but degen- 
erates to K y = K z = 0, where k p = (K x , K y , K z ), 
k q = (K — K x , —K y , —K z ). There are three cases. 

1. For K x < < K, |kp| = -K x , \k q \ = K + K x , 

Sp S, Sq S. 

2. For < K x < K, \k p \ = K x , \k q \ = K — K x 
s p = s q = s. 

3. For < K < K x , \k p \ = K x , \k q \ = K x — K, 

S p 5, S q S. 



Let us expand the Hamiltonian (2.11) (expressed in 
terms of b, b*) in power series: 



Tt = Ti-o + "Hit 



(2.17) 



Here H.q is quadratic in b and b* , giving the Hamiltonian 
of non-interacting waves: 



Wo = y ckb{k)b*{k)dk, 



(2.18) 



with linear dispersion relation tt>o(k) = ck. In the Hamil- 
tonian of interaction Tl m t we take into account only three- 
wave processes: 



G 



Hint = ~ / (^(k,k 1 ,k 2 )6^ 2 6 3 + c.c.) 

X(5(ki -k 2 - k 3 )dkidk 2 dk 3 . (2.19) 

We neglected here 3 processes (processes described 
by 6*6263 and 616263 terms), because they are nonreso- 
nant. It means, that if we take into account <-> 3 term, 
it is not going to change our final results, thus we can ne- 
glect it from the very beginning. We also neglected con- 
tributions from 4- wave and higher terms, because three- 
wave interaction is the dominant one. 

The coupling coefficient of the 3-wave interaction a 
given by 



V(ko,ki,k 2 ) 



(2.20) 



(3<7 + COS 0oi + COS O 2 + COS 012 ) 



where g is some dimensionless constant of the order of 
unity and 9ij is the angle between k^ and k,. Since we 
have almost linear dispersion relation, only almost par- 
allel wavevectors can interact, therefore cos 9 a with the 
high accuracy can be replaced by 1 and ( [2.20 ) reduces to 



y(k ,ki,k 2 ) 



' ckk\k2 



3(9 + 1), 



(2.21) 



3. Canonical Equation of Motion 



The Hamiltonian equations of motion ( 2.1 2| ) for the 
complex canonical variables 6 6* have standard form [EJ 



,db(k,t) _ SH 



dt 



(2.22) 



<56* (k, t) 

For the acoustic Hamiltonian ( 2 . 1 7 - pTT9| ) . this equation 
takes the form 



4t ck . 



6(k,t) 



= \J V(k,q,p)6(q)6(p),5(k-q-p)^ P 
+ j y*(k,q,p)6(q)*6(p)£(k + q-p) 



(2.23) 



d 3 qd 3 p 
(2tt) 3 ' 

It is sometimes convenient to concentrate attention on 
steady state turbulence, which is convenient to describe 
in the k, u> -representation. After performing a time 
Fourier transform, one has instead of (2.2S), 

c*]&(k,w) = \ J F(k,ki,k 2 ) 

dkidLUidk.2duj2 



x 6i&2<5(k — ki — lt2)S(u> — u>i — L02) 



(27T) 



- J V*(k, ki, k 2 )b* 1 b 2 S{k + ki - k 2 )5(u wa) 

(2.24) 



dkidujidk2duj2 



(2tt) 



Hereafter we will refer to this as the basic equation of mo- 
tion for the acoustic turbulence normal variables bf., b* k 
and use it for a statistical description of acoustic turbu- 
lence. 



C. Relations between Wave Amplitudes 
a + (k), a~(k) with Normal Variables of Acoustic 
Turbulence b(k), fo*(k). 



Comparing Eqs. (2.4) and (2.5) we get 



8p(k,t) = poe|a + (k, t) exp[iuj(k)t] 

+ a- (k, t) exp[-Kj(k)i] j (2tt) 3/2 

$(M) = ^{o + (k,t)exp[t«(k)t] 



(2.25) 



a~ (k, t) exp[-iu{k)t] } (2tt) 3/2 (2.26) 



Here $ is velocity potential: v = V$. This gives 
exp[— iLo(k)t] 



a + (k,t) 
a~(k, t) 



2e(27r) 3 / 2 

exp[iw(k)t] 
2e(27r) 3 /2 



« P (k,t) 



•i*(M)^ 
cr 



(2.27) 



Note, that a + and a is dimensionless variables. 

Now we can easily express a + (k), a'k) in terms of 
6(k), 6*(k) and thereby relate the two alternative ap- 
proaches presented in this paper, 



a + CM) - - W — (27r)- 3 / 2 exp[-^(k)t]6*(-k), (2.28) 
e V 2cp 



a~(k,*) = - \ — (27r)- 3/2 exp[ia;(k)t]6(k) . (2.29) 
e y 2cp 



To check, that the two approac hes are consistent, we 
rewrite the equation of motion (2.6) for a s k neglecting 
e 2 (four-wave interaction) terms: 



da s (k,t) 
dt 



= ^E/ dk p dk q L s ^Xa^(k pi t)a^(k q ,t) 



x 5(k p + k q — k) exp{i[s p cj(kp) + s g w(k 9 ) — sw(k)]t} 



(2.30) 



Now we substitute Eqs. fl2.28| ) and ( p3s| ) into fl2.30| ) and 
obtain 
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dt 



iw(k) 



b(k,t) 



dpdq* 



' kpqc 
4n 3 p 



teroo'O" ^k,-k',-k"9k',-k'?k",-k" 



(// — 2) + COS 0k, p + COS 0k, q + cos ( 



p,q 



(2.31) 



exp[z(s p u/ + s q ui h 



>)t], 



k + k' + k" = 



(3.6) 



\5(k + p + q)b*b* + 25 k+p _ q 6;6 q + 6Ck - p - q)& p fe q l where the symbol P 00 ' (Poo'o" ) means that 
LV ' p q P q p q V all replacements -> 0', 0' -> ( 0', 0' -> 0", 



we sum over 



Now one can see that eq uatio n (2.31) looks exactly as 
( 2.22 ) with Hamiltonian ( [2. IS ) and with coupling coeffi- 
cient (2.20). Thus one conclude that the two approaches 
are equivalent and the choice between them is the ques- 
tion ol taste. 



3. LONG-TIME ANALYSIS OF STATISTICAL 
BEHAVIOR 

The analysis proceeds by first forming the hierarchy of 
equations for the spectral cumulants (correlation func- 
tions of the wave amplitudes) defined as follows. The 
mean is zero. 



a s (k)a s (k')) = S(k + k')q ss (k, k'), (3.1) 

a s (k)a s '(k')a s "(k") 

5(k + k' + k")q ss ' s "(k, k', k"), (3.2) 
^a s (k)a s '(k')a s "(k")a s "'(k"')) 

<5(k + k' + k" + k"')q ss ' s " s '" (k, k', k", k'") 
+<5(k + k')<5(k" + k"')q ss ' (k, k')q s " s "' (k", k'") 
+8{k + k")6(k' + k"')q ss " (k, k")q s ' s "' (k', k'") 

5(k + k"')<5(k' +k")g ss '"(k,k'")/ s "(k',k"), (3.3) 



where (...) denotes average and the presence of the delta 
function is a direct reflection of spatial homogeneity. In- 
deed the property of spatial homogeneity affords one a 
way of defining averages, which does not depend on the 
presence of a joint distribution. We can define the av- 
erage (^(x)rj(x + r)) as simply an average over the base 
coordinate, namely 



(r?(x)r?(x + r)) 



1 



{2Lf 



?7(x)?7(x + r)dx . (3.4) 



-L 



To derive the main results of this paper, it is sufficient 
to write the equations for the second and third order cu- 
mulants. They are 



dt 



e-xp[i(s p uj p + s q uj q - sui k )t]8(k - p q), 
k + k' = ; 



0, -> 0",0' -> 0, 0" -> 0'). 

The total energy of the system per unite volume can 
be written as 

lirn(ip u j (x)« i (x + r) 

r— >0 Z 

2 2 
C On PC 

+ — + r) + - 2)r ? (x)? ? (x + r ) 



limV 



2 1 



(l- Sl s 2 )g SlS2 (k)e jkr dk 



lim / p c z e'(q + -(k) + q- + (k))e lKr dk 

r— >0 



= / 2pocVg+-(k)dk 



(3.7) 



since q + (k) = q + (— k). The spectral energy is there- 
fore 2pc 2 e 2 g H (k). For convenience we denote q^ (k) as 
e(k). 

To leading order in e, q ss (k, k') and q ss s (k, k',k") 
(which we may call qQ S (k, k') and q^ s s (k, k',k")) are 
independent of time. Anticipating, however that certain 
parts of the higher order iterates in their asymptotic ex- 
pansions may become unbounded, we will allow both 
qQ S (k, k') and qg s s (k, k',k") to be slowly varying in 
time 



dqf (k,k>) 
dt 

dg s s ' s " (k,k',k") 
dt 



c 2 tjss 1 



2 jpss' s" 



(3.8) 



and we will choose F2 and F3 to remove those terms with 
unbounded growth from the later iteration. We will find 
that for s' — — s, F^ s is given by the right-hand side of 
acoustic KE: 



TTiSS S S 

b 2 - % 



xq Sp ' Sp (k p )A(s p uj p + s q u q - SLu)S(k p + kq - k)dk p dkq 



lim / SqSp f (L s ' sp,Sq V 



:q° s (k q , -k q )g s *' s " (k p ) A(s p aj p + s q uj q 



(3.5) 



S(k p + kq — k)dkpdk q 



eP 0Q , „ J dkpdk q L^X^%X 6 ^ " P " 1) ^ that ^ ^ F "''" 



dt 

X exp[i(s p LU p + SqUJg - SUJ k )t] + 



have the form 



iC(k,k')(O k + n k ,) 



(3.9) 
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and 



iq^ ss (k,k',k")(0£ 



(3.10) 



respectively. It is clear that Cl s k can be interpreted as a 
complex frequency modification. Its exact expression is 
given by 



nt = -4i lim V 



(j-^S,S p ,S q \ 



(3.11) 



x <j Sp ' Sp (kp)A(s p ujp + s q u q - suj)S(k p + k q - k)dk p dk q 
and, whe n calc ulated out, is preci sely equal to s(ui 



c\k\)e 2 in Jl 171) ■ Note that in ( |3 . , i = Te 2 and T is fi- 



nite. The ln(l/e ) cocfhcicnt comes from the term hit or 
ln(T/e 2 ) = InT + ln(l/e 2 ) in the asymptotic expansion. 
For finite T, the dominant part is ln(l/e 2 ). 

The perturbations method has the advantage that it is 
relatively simple to execute. However, there is no a priori 
guarantee that terms appearing later in the formal series 
cannot have time dependencies which mean they affect 
the leading approximations on time scales comparable to 
or less than e -2 (e.g. a term e 4 t 3 should be accounted 
for before the term e 2 t). To check this, one must have 
a systematic approach for exploring all orders in the for- 
mal perturbation series and removing (renormalizing) in 
groups those resonances which make their cumulative ef- 
fects at time scales e" Ar (ln(i)- M ), N, M = 1,2,3,.... 
The diagram approach, which requires some familiarity 
to execute, is designed to do this and, both for complete- 
ness and the fact that we will have to proceed beyond 
the one-loop approximation to resolve the questions of 
the angular redistribution of spectral energy, we include 
it here. 



We will be interested also in the double correlation func- 
tion n (k, u>) of the acoustic field b, b* 

(27r) 4 n(k, u)S(k - k')5(u> - u') = (6(k, Lo)b*(k', to')) . 

(4.3) 

The simultaneous double correlator of the acoustic field 
n(k) is determined by 



(27r) 3 n(k)^(k-k') = (&(k,i)6*(k',i)) 



(4.4) 



This is related to the different-time correlators in the u> 
representation n(k, lu) as follows: 



n(k) = / n(k, lu) 



2^ 



(4.5) 



The Green's and correla tion functions together with 
the bare vertex V(k, q, p) ( 2.2C ) are the basic objects of 
diagrammatic perturbation approach which we are going 
to use (see Fig la). 



FIG. 1. Panel (a): Basic objects of diagrammatic pertu- 
bation approach. Panel (b): First terms in the expansion of 
mass operator E(k, ui). 



FIG. 2. Diagrams (a) from Fig.l with specified directions 
of arrows. 



4. DIAGRAMMATIC APPROACH TO 
ACOUSTIC TURBULENCE 

A. Objects of Diagrammatic Technique 



Let us define the "bare" Green's function of Eq. (2.24) 



a-. 



G (k) 



f 



lu — ck + iO 



(4.1) 



One may see from (2.24) that this function describes the 
response of the system of noninteracting acoustic waves 
on some external force. We added in the denominator the 
term +i0 by requirement of causality. We remark that 
causality (the arrow of time) is introduces in the pertur- 
bation approach by the limit t — * oo and the fact that 
sgni appears in (1.7). Next we introduce the "dressed" 



Green-function which is the response of interacting wave 
systems on this force: 



(27r) 4 G(k, oj)S(k - k')6{uj - lu') 



I Sb(k,uj) 

W) 



W(k' 



(4.2) 



B. The Dyson- Wyld equations 

In the diagrammatic series for the Green's function one 
may perform the partial Dyson's summation over one- 
particle irreducible diagrams. This results in the Dyson 
equation for the Green's functions: 



G(k,«) 



f 



lj - (jjoQt) -MO - £(k,w) 



(4.6) 



where the "mass operator" E(k, to) gives the nonlinear 
correction t o the complex frequency Wo(k) + i0 due to the 
interaction ( 2.19 ). This is an infi nite s eries with respect 
to the ba re am plitude V(k, q, p) (2.2C), dressed Green's 
function (4.2) and double correlation function n(k, to) 
(4.3). All of the contributions of the second and fourth 
order in V are shown on Fig 1(b). 

We have not specified the direction of arrows on Fig 
1(b); each diagram should be interpreted as a sum of 
diagrams with all possible directions of arrows compati- 
ble with vortex V(k, q, p), describing the three- wave pro- 
cesses 1 <-> 2. For example, diagram (a) on Fig 1(b) cor- 
responds to three diagrams shown on Fig 2. The diagram 
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(a4) on Fig 2 describes the nonresonant process 0^3 
which is not essential for our consideration. 

With the help of the similar Dyson's summing of 
one-particle irreducible diagrams, one can derive Wyld's 
equation for n(k, ui): 



n(k, w) = |G(k, u)\ 2 [£>(k, ui) + $(k, cj)] 



(4.7) 



Here D(k, w) is the correlation function of white noise, 
(27r) 4 D(k, u>)5(k - k')5(u> - J) = (f(kw)f*(k'uj')) , 

(4.8) 

and the mass operator <j>(k, a;) describes the nonlinear 
corrections to D(k, ui). This is an infinite series with re- 
spect to the same objects G(k, u>),n(k, ui) and V(k, q, p). 
All diagrams of the second and fourth order are shown 
on Fig 3(a). 

We also have not specified arrow directions in the di- 
agrams for S(k, ui) and 3>(k, ui). In complete analogy 
with diagrams for G(k, ui) one diagram on Fig 3a corre- 
sponds to two diagrams (al) and (a2) on Fig 3b. All the 
rest diagrams for $(k, ui) reproduces in the same way - 
one chooses all possible directions of arrows and discards 
those which incompatible with definition of vertex V (see 
Fig la). 



FIG. 3. First terms in the diadrammatic pertubation ex- 
pansion for mass operator ^(k^). 



w(k) =w (k) +Re£(k,w* 
7(k) = -ImS(k,w*) . 



(4.12) 
(4.13) 



Now we have to decide how to choose "in the best 
way" . The sim plest way is to put uj* = cJo(k) = ck, as it 
was stated in ( l.lOj ). As a next step we can take "more 
accurate" expression lu* = u>(k), i.e. to take into account 
the real part of correction to Wo(k). But later we will see, 
that better choice is 



cj* = ui(k) + *7(k) 



(4.14) 



which is consistent with the position of the pole of 
G*(k, ui). We will show that this choice is self consis- 
tent while deriving the balance equation in section 5.3. 



2. The double correlation function 

The same type of approximation may be performed 
for t he c orrelation function. Namely in the Wyld equa- 
tion d4.7| ) one may replace G(k, ui) by G(k, ui) and to ne- 
glect the ui dependence of <I>(k,w) by putting $(k, ui) — > 
$(fe) = $(k,w*), or 



n(k,w) = |G(k,w)| 2 D(k) + $(k) 



(4.15) 



We will call this one-pole approximation for the correla- 
tion function. 



D. One- loop approximation 



C. One-pole approximation 

1. The Green's function 

We have assumed from the beginning, that the wave 
amplitude is small. Therefore, 



E(k,w) < ui (k) . 



(4.9) 



As a result the Green's function has a sharp peak in the 
vicinity of ui = ck and one may (as a first step in the 
analysis) neglect the w-dependence of E(k, u>) and put 



S(k,w) ~ S(k,w ~ ck) 



(4.10) 



The validity of this assumption will be che cked later. Un- 
der this assumption the Green's function ( |4.2| ) has a sim- 
ple one-pole structure: 



<5(k,w) 



ui — w(k) + «7(k) ' 



(4.11) 



where 



Let us begin our treatment with the simple one-loop 
(or direct interaction) approximation for mass operators 
£ and <F This approximation corresponds to taki ng int o 
account just the second order (in bare vertex V ( [2.20 )) 
diagrams for the mass operators E and Two loop ap- 
proximation will be considered in Appendix C. We will 
estimate two-loops diagrams and we will show, that some 
of them gives the same order contribution to 7^ as one- 
loop diagrams. Therefore, one loop approximation is an 
uncontrolled approximation, but we believe, that it gives 
qualitatively correct results. Note that these diagrams 
include the dressed Green's function in contrast to the 
approximation of kinetic equation which is nothing but 
one-loop approximation with the bare Green's function 
inside. We will see later that this difference is very im- 
portant in particular case of acoustic turbulence. The 
KE for waves with linear dispersion law forbids the an- 
gular evolution of energy because conservation laws of 
energy and momentum allow interaction only for waves 
with parallel wave vectors. In the one-loop approxima- 
tion with dressed Green's function, the conservation laws 
ui(k) ± w(ki) = ui(k ± ki) are satisfied with some accu- 
racy [of the order of 7(k)]. As a result, there exists a 
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cone of allowed angles between k and ki in which inter- 
actions are allowed. Therefore one has to expect some 
angle evolution of wav e pack ages w ithin this approxima- 
tion. Combining ( 4.11 ) with ( |4.15| ) one has the following 
expression: 



n(k, lu) 



2 7 (k)n(k) 



[u;-Lu{k)} 2 + 1 2 (k) 



(4.16) 



1. Calculations o/E(k,o;) 

In the one- loop approximation expression for E(k, lu) 
has the form 

E(k, lu) = E„i(k, lu) + E a2 (k, lu) + E a3 (k, u) , (4.17) 

where Ej(k, cj) is given by (|Al]-|A3|). Our goal here is to 
analyze these expressions in one-pole approximation, by 
substituting in it "one-pole" ?i(k, lu) and G(k, lu) from 
(4.11) and ( 4.16| ). In the resulting expression one can 
perform the integration over lu analytically. The result is 



E(k,w) 



d 3 k 1 d 3 k 2 



(2tt) 3 

|V(k 2 ,k,ki)| 2 J(k + ki-k 2 )[n(ki 



(4.18) 



V uj + w(ki) - w(k 2 ) + i(7i + 72) 

|y(k ,k 1; k 2 )| 2 J(k-k 1 ~k 2 )n(k 2 ) \ 
cj - w(ki) - w(k 2 ) + 1(71 + 72) I 

Next we introduce E(k) — E(k, lu*) , with lu* given by 



(4.14) and consider (4. IS) in the limit of small 7, which 
allows us to perform analytically integrations over per- 
pendicular components of wavevectors. The result for 
the damping frequency "/(k) may be represented in the 
following form (for details see Appendix B): 



7 (fc) 



A 2 k 2 
Aire 



2 1,2 



l/L 



n(q)q 2 dq 



A 2 k 



N(Q) 



(4.19) 



We introduced here cut-off for small k at l/L, where L 
is the size of the box. We also introduced "the density of 
the number of particles" N(Q) in the solid angle accord- 
ing to 



N(Cl) = I k 2 n{k)dk , 
such that the total number of particles 

.V / N(n)dn . 



(4.20) 



(4.21) 



After substituting A from (Bll), one has the following 
estimate for 7(k): 



7 (k) ~ k 2 N(n)/ Po , 



(4.22) 



Consider now E'(k) 
that 

A 2 



E(k) = 



47T 2 C 

x[(fc 2 
A 2 k 



ReE(k). It follows from flBl^ ) 

(4.23) 



dq I dyq n(q) 


2 



2kq + q 2 ) 
dq 



(k 2 
ydy 



r 2 

1 fel2 



r 2 

1 fcl2 



2kq + q 2 )] 
[cq 3 n{q)] , 
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where r^i 2 = j(k) + 7(^1) + 7(fc 2 ) is the "triad inter- 
action" frequency. One may evaluate the integral with 
respect to y as 

L ( 5 ) =m ^~l n ^L (4.24) 



kkq 



Ti (k) ' 



After substituting 7(k) from ( 4.22 ), one has 

L(q) oclnp /giV(n) (4.25) 

The main contribution to the integral (4.23) over q comes 
from the infrared region q ~ l/L. It gives the estimate, 



S'(k) = 4^Li?(fi) 



(4.26) 



where we have defined the density of the wave energy in 
solid angle as 



E(Sl) = J Lu a (k)n(k)k 2 dk . 

This value relates to iV(fi) as follows: 
E{n) 



(4.27) 



(4.28) 



Equation ( 4.26 ) together with the expression ( Bll ) for 
A may be written as 

E'(k) ~ cfcelnl/e , (4.29) 

where 

e ~ E(n)/p Q c 2 , (4.30) 

is the dimensionless parameter of nonlincarity, the ra- 
tio of energy of acoustic turbulence and the density of 
thermal energy of media poc 2 ~ hT, where n is the con- 
centration of atoms. 



Equation (4.22) for 7(k) may be written in a similar 
form 



7(k) ~ ck(kL)e 



One can see that 



7 (k) 



kL 



oc 



(4.31) 



(4.32) 



E'(k) lnl/e ' 

It means that for large enough inertial interval 

7 (k) > E'(k) (4.33) 

and one may neglect the nonlinear corrections E'(k) to 
the frequency with respect to the damping of the waves 
7(k). That shows that our above calculations of E(k) is 
self-consistent. Later we also will take into account only 
damping 7(k) in the expressions for the Green's functions 
taking u){k) = LUo(k) = ck. 



11 



2. Calculations of &(k,u) . 

In the one-l oop approximation expression for <I>(k,w) 
has the form (A4). Afte r sub stitution of n(k, ui) in the 
one pole approximation ( 4.16 ) one may perform analyti- 
cally integration over frequencies: 



$(k,w) 



d 3 k 1 d 3 k 2 ,,,,,, 
(2tt)3 " (kl)n(k2) 



(4.34) 



\V(k, kj, k 2 )| 2 ( 7 (k 1 ) + 7 (k 2 )),5(k - ki - k 2 ) 
[w - w(kO - ^(k 2 )] 2 + [ 7 (k!) + 7 (k 2 )] 2 

\V(k 2 , ki, k)| 2 ( 7 (k!) + 7 (k 2 ))(S(k + ki - k 2 
[to + w(k!) - w(k 2 )] 2 + [ 7 (kx) + 7 (k 2 )] 2 



We will analyze this expression in the next section. 



3. Balance Equation 



Consider the Dyson- Wyld equations ( i£ ) and ( |4.7|) in 
the inertial interval, where one can neglect 7 o(k) in com- 
parison with ImS(k, ui) and D(k) in comparison with 
$(k,w): 



G(k,w) 



1 



w — cjo(k) — S(k, w) 
T»(k,w) = |G(k,w)| 2 $(k,w) . 



It follows from (4.6) that 

ImG(k,w) = |G(k,w)| 2 ImE(k,w) 



(4.35) 
(4.36) 

(4.37) 



By comparing ( 4.36| ) with ( 4.37 ) one may see that the 
following combination 

L(k, lu) ee $(k, w)ImG(k, w) - n(k, w)Im£(k, ui) (4.38) 

is equal to zero. In particular 



L(k) = | L(k,c)^=0 



(4.39) 



Together with ( 4.3S ) it gives 

Im / ^ [$(k,w)G(k,u;) - n(k, w)E(k, w)] =0 . (4.40) 

Let us co mput e now the first term i n (|4.40| ) . By substitu- 
tion Eq. ( 4.34 ) for $(k, ui) and Eq. ( |4 . 1 1| ) and integration 
over u) one has 



G(k,w)$(k,w) 



^kirfk 2 n , n , 
— n(ki)n(k 2 ) 



^ _ v _,_, - v _,_, j ^ 

1 l^kLkajl^k-ki-^) 

2 w (k) - w (ki) - w (k 2 ) - ir fe i2 

^(ka.ki.kjp^k + k! -k 2 ) 
wo(k) + w (ki) - w (k 2 ) - iT kl2 



(4.41) 



Next we will perform integration over u> in ( 4.40| ). Re- 
member S(k, a;) is analytical function in the upper half 
plane of u> while n(k, ui) has one pole there. Therefore 



Im / ^n(k,u>)E(k,uj) = n(k)Im£(k,w*) (4.42) 

J 27T 



where is given by ( 1.15 ). This is the justification of 
our choice lu* . 

Now let us put everything together to obtain 



= L(k) = | ^^r^j^k-kx-k,) (4.43) 

k 1 |y(k,k!,k 2 )| 2 [n(ki)Ti(k 2 ) - n(k)[n(k!) + n(k 2 )]] 
X 2 (w (k) - wo(ki) - coo(k 2 )) 2 + T 2 12 

+ ^(k + k!-k 2 ) 

|V(k a ,ki,k)| 2 [n(k a )[n(ki) + n(k)] - n(k)n(ki)] 1 



( Wo (k) + ^ (k 1 )- Wo (k 2 )) 2 + r 2 
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This is the main result of the diagrammatic approach: 
the balance equation for stationary in time acoustic tur- 
bulence. In nonstationary case one can get similarly the 
generalized kinetic equation in the form 



dn(k,t) 
dt 



L(k,t) 



(4.44) 



where L(k, t) is given by Eq. ( 4.43| ) with correlator de- 
pending on time n(kj) — > n(kj,t). In the limit 7 (k) — > 
this expression turns into well known (cf. Q) collision 
integral for 3-wave kinetic equation 



dki<ik 2 

3~ 



St\n(k,t)\ = lim L(k) = 2vr f ^ 

1 V n r fcl2 ^o v ' J (2tt 

x i,5(k-k 1 -k 2 )|y(k,k 1 ,k 2 )| 2 

{n(k x )n(k 2 ) - n(k)[n(ki) + n(k 2 )]} 
x S[ui (k) - wo(ki) - w (k 2 )] 
+ ^(k + k x -k 2 )|T/(k 2 ,kx,k)| 2 

{n(k 2 )[n(kx) + n(k)] - n(k)n(ki)} 
x tf[wo(k) + wo(ki) - woM • (4-45) 



We see that the generalized kinetic equation differs from 
the well known collision term in the three wave kinetic 
equation by replacing 5-functions on the corresponding 
Lorenz function with the width of r,ti 2 -triad interaction 
frequency. 



12 



5. CONCLUSION 



In the present paper we have begun to develop a consis- 
tent statistical description of acoustic turbulence based 
both on the long-time asymptotic analyses (Section 3) 
and on the perturbation diagrammatic approach (Sec- 
tion 4). The first approach is more straightforward. The 
diagrammatic approach provides a systematic way of an- 
alyzing higher order terms in the perturbation theory. 

Our main result is that the nonlinear corrections to 
the frequency is much smaller than the nonlinear damp- 
ing of the waves. We find also the balance equation (4.43) 



which generalizes the simple kinetic equation for acoustic 
waves. One can show that the balance equation (4.42) 
has the same isotropic solution (Zakharov-Sagdeev spec- 
trum) as the kinetic equation. However the kinetic equa- 
tion for acoustic turbulence does not describe the angle 
evolution of turbulence: any arbitrary angle distribution 
i s the solution of KE. In contrast, our balance equation 
(4.43) contains terms which describe an angular redistri- 
bution of the energy because of the non-zero value of the 
interaction cone, which is proportional to Tku- But we 
have yet to show that this expression contains all such 
terms to this order. 

One may imagine three very different ways of the an- 
gle evolution of anisotropic acoustic turbulence. The first 
one is a tendency to form very narrow beams with the 
characteristic width of about interaction angle. The sec- 
ond one is an approach to isotropy downstream to the 
large wave vectors. The last possibility is to form a beam 
with a characteristic width of about unity, exactly like it 
happens in the turbulence of waves with weak dispersion 
|l7| . Another important question is: does the spectra 
of acoustic turbulence depend on the features of pump- 
ing or they are universal (independent of details of en- 
ergy influx)? We intend to answer these questions (in 
the framework approximations we made in that paper) 
in our next project. 

It is an exciting challenge to try to go beyond the ap- 
proximations made here in order to understand whether 
the scaling index of the interaction vertex in the system 
of acoustic waves in two- and three-dimensional media 
must be renormalized or not. 
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APPENDIX A: RULES FOR WRITING AND 
READING OF DIAGRAMS FOR MASS 
OPERATORS 

Here we state without proof the set of rules for writing 
down diagrammatic series: 

1. In order to write down all diagrams for E and $ 
of 2n order in vertices, one should draw 2n vertices 
and connect them with each other by lines n and 
G in all possible ways. Two ends must be left free. 
If both ends are straight, we shall get a diagram 
for ^(k, a;); if one of them is wavy, this will be a 
diagram for E(k, u>). 

2. The diagrams for $ and E containing closed loops 
in GF are absent. This follows from the fact that 
the Wyld's DT appears from glued trees. 

3. There is no mass operator with two wavy ends in 
DT. 

4. In the diagrams for $ (for E) one can pass from 
every vertex along the G lines to the entrance and 
to the exit in a single way. 

5. In every diagram for E there is a single root link- 
ing the entrance and exit along the G lines - the 
backbone of the diagram. The rest G lines of the 
diagrams may be called the rips. 

6. The diagrams for <E> contain the basic cross section 
in which they may be cut in a single way into two 
parts only at lines n(k, lu). 

7. Every V vertex is entered by one arrow and excited 
by two. The V* vertex is entered by two arrows and 
exited by one. 

One can show (see |l2]] ) that rules (3-7) follows from 

The rules of reading diagrams are the follows: 

1. Write down product of DT objects (double correla- 
tor, Green function or vertex) (with corresponding 
arguments) corresponding to each element of the 
diagram. 

2. Write down delta- functions in 4-momenta for 2n — 1 
vertices in such a way, that the sum of entering 4- 
momenta is equal to the sum of exciting ones. One 
of vertices (for example the one corresponding to 
the end of the diagram) does not contain the delta 
function. 

3. Perform integration along all internal lines of dia- 
gram : di — dki/(2Tr) d du>i/ (2ir) where d is space 
dimension. 

4. Then you have to multiply diagram by (2ir)( d+1 \ 

5. To multiply diagram by l/p where p is the number 
of elements in its symmetry group. 
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For example diagrams (al), (a2) and (a3) correspond to 
the following analytical expressions: q 



S a i(k,w) = J 



d 3 k\d 3 k 2 duj\duj2 



(2tt) 3 2tt 
x S{uj + uj 1 - w 2 )|^(k 2 ,k,k 1 )| 2 G'2ni, 



5(k + ki-k 2 ) 

(Al) 



E o2 (k,o;) = y -— 



d 3 k 1 d 3 k 2 dw 1 dw 2 , , , N 
(2^ ^(k + fc-k,) 



x + - w 2 )|y(k 2 ,k,ki)| 2 G 1 n 2 , 
d 3 k\d 3 k 2 du>idu>2 



(A2) 



/ d i k 1 d i k 2 dw Y du 2 xn , , , 
S 3(k,a;) = y -— ^ __ 5(k _ kl _ k2 ) 



(2tt) 3 2tt 
x £(w-Wi - w 2 )|y(k,ki,k2)| 2 Gin 2 , 



(A3) 



We defined here the following shorthand notation, Gj = 
G(ky, u>j); n, = 7ii(ki,<x>i) In the same way one can find 
analytical expressions for (f> a (k, cj): 

f d 3 kxd 3 k 2 du>iduj 2 I 2 

$(k,^) = / -7^3 ^— [ol^( k > k i: k 2)r 



(2tt) 3 2tt l 2' 
x riin 2 (5(k — ki — h 2 )S(uj — lj\ — lj 2 ) (A4) 
+ |F(k2,ki,k)| 2 ni^2^(k + ki - k 2 )S(uj + ui - w 2 )] ■ 

Analytical expressions for s the 4-th order diagrams (two- 
loop diagrams) will be shown in Appendix C. 



APPENDIX B: CALCULATION OF £(K,w) 
-DETAILS 



Let us start from ( LIE ) and introduce E(k) = E(k, w*) 
with given by (4.14): 



E(k) 



d 3 kid 3 k 2 
(2tt) 3 



(Bl) 



/ |y(k 2 ,k,k 1 )| 2 J(k + k 1 -k 2 )[n(k 1 )-Ti(k 2 )] 
y w(k) + w(ki) - w(k 2 ) + ir fe i2 

|V(k , kx, k 2 )| 2 5(k k x k 2 )n(k 2 ) \ 



where 



w(k) - w(ki) - w(k 2 ) + ir fc i 2 



r fe i 2 = 7 (k)+7(ki) + 7(k 2 ) 



(B2) 



is "triad-interaction" frequency and 1/Tfci2 is triad inter 



action time. One can consider (Bl 



i) as an integral 

equation for the damping of wave 7(k) = — ImE(k) and 
for the frequency w(k) = o->o(k) + ReE(k). 

First we consider these equations in the limit of weak 
interaction where , T — > 0, and the main contribution to 
the first term in (Bl) comes from the region where 



w(k) + w(ki) = w(k 2 ), 



k + ki 



(B3) 



These are conservation laws for 3-wave confluence pro- 
cesses + 1 — ► 2. The main contribution for the second 



term in ( B 1 ) comes from the region 



w(k) = w(ki) + w(k 2 ), 



ki 



(B4) 



These are conservation laws for decays proces ses — ► 
1 + 2. For weak interaction one may replace in (4.2) and 
( p~3| ) w(k) on uj Q (k) = ck. Than it follows from ( p~2] ) 
and ( |4.3| ) that ki || k 2 || k with ki,k 2 directed along 
k. This fact makes it natural to introduce in integrals 
(Bl) new variables: scale positive variable q > and 
two-dimensional vector k such that 



ki = qh/k + k , k _L k 



In the first term of (Bl) 

k 2 = (k + g)k/fc + k, 
In the second term 

k 2 = (fc — g)k/fc — k, 



< q 



0<q<k 



(B5) 



(B6) 



(B7) 



For k <C k the denominators in integrals (Bl) strongly 
depend on k. Indeed: 



cd (k) +w(ki) - w(|k + ki|) = cfc 7 



(B8) 



2g(fc + <z) ' 

w (k)-o;(k 1 )-a;(|k-k 1 |) = -cfe- 7 ^ -. (B9) 

2g(fc - q) 

This allows to neglect k dependence of interaction 
V(k, q, p) and correlation n(ki) in numerator of (Bl) for 
estimation. The result is 



A 2 k 



dn x 



dq 



dq- 



g{k + q) [n(q) - n(k + q)] 
ckn 2 /[2q(k + q)} + jr fc i 2 

q(k-q)n(q) 



-ckK 2 /[2q{k~q)]+ir k 
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where 



A = 3(«?+l) v /c/47r 3 p , 



(BIO) 



(Bll) 



is a factor in ( |2.20| ) so that for parallel or almost paral- 
lel wavevectors I^(k, q, p) = A^/kqp. After changing of 
variables this integral becomes to be more transparent: 



E(k) = 



A 2 

47T 2 C 



dq 



- njk + q)] 
y + 7r fe i 2 



dyq 2 (k + q) 2 
dy 



(B12) 



fc q 2 {k-q) 2 n(q) 
dq 

y - «r fe i 2 



One may estimate 7/ max — ck 2 /2q from the fact that our 
expressions were obtained by expanding in n/k, therefore 
should be at least k < k. 

Now let us consider the imaginary and real part of 
E separately. It is convenient to begin with 7(k) = 
-ImE(k): 
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7 (k) 



A 2 

AlT 2 C 



r OO pOC 

dy dqq 2 (k + q) 2 r kl2 (B13) 
Jo Jo 



[n(q) - n(k + q)] 
ii 2 4- r 2 

y + 1 fci2 



dq 



q 2 (k-q) 2 n(q)T kl 

a 2 + r 2 
y + 1 fci2 



Here we changed the upper limit of integration: y max — ► 
oo because the main contribution to the integral comes 
from the area y ~ T <C ck. After trivial integration with 
respect of y one has: 



A 2 



q 2 (k + q) 2 [n(q] - n(k + q)] dq 



q 2 (k - q) 2 n(q)dq 



(B14) 



This expression for 7(k) corresponds to that given by 
the kinetic equation jlj for waves. For further progress 
it is necessary to do some assumption about n{q). Let 
us assume that n{q) vanishes with growing of q faster 
than 1/q 4 . 1 For such spectra the main contribution to 
the integral comes from small q <C k. In t his c ase contri- 
butions from first and second integrals in (B14) coincides 
and may be represented in the form: 



7 (fc) 



A 2 k 2 

AlTC 



1/L 



n(q)q 2 dq ■ 



A 2 ^ 
Aire 



N(Q) 



(B15) 



In this c ase contributions from first and second integrals 
in (B14) coincides and may be represented in the form: 



7 (fc) 



A 2 k 2 
Aire 



n(q)q 2 dq • 



1/L 



A 2 ^ 
Aire 



iV(fi) 



(B16) 



APPENDIX C: ESTIMATION OF THE 
TWO-LOOP DIAGRAMS. 



Let us write down analytical expression which corre- 
spond to one of the diagrams (b) in Fig. 1(b) 



E 6 (k,w)= J 



dkidk2<Ai>iG?a->2 
(2nf 



V a V b V c V d n(ki, wi)n(fc 2 , u 2 ) 



v a = y(k 1 + k 2 + k,k,k 1 + k 2 ), 

V r 6 = V(k 1 + ka+k I ki > k + k 3 ) J 

y c = y(k 2 + k,k 2 ,k), 

V r d = V(k 1 +k a ,k 1 ,k a ) 



(C2) 
(C3) 
(C4) 
(C5) 



We just followed the rules of DT and integrated over all 
delta-functions. From now, the an alyses will be parallel 
t o tha t of appendix B. Let us use ( 4.16 ) for n(k,u>) and 
(4.11) for G(k, oS). Now we can easily perform integration 
over u)\ and uj 2 . Now, as it was done in appendix B, intro- 
duce £&(k) = S(,(k, w*). Since all interacting wavevec- 
tors are almost parallel, we introduce two-dimensional 
vectors K\ and k 2 such that 



ki = gik/fe + K\ , Kt A- k 
k 2 = q 2 k/k + k 2 , k 2 A- k 



(C6) 
(C7) 
(C8) 



We use V(k, q, p) = A\fkq~p. Si nce m <C k, we can ex- 
pand resonance denominators in (CI) with respect to /t,. 
The integrals will be dominated by regions, where qi <C k. 
By putting everything together, one gets 



-K 2 dq x dq 2 dK 2 dK 2 4 3 . _ 

—A k [qi + q 2 )q\q 2 n qi n q 



2tt 6 



2(91+92) 



91,92,91+92 



)( 



2 qi 



<l2 



ilk) 



2q 2 



ilk) 



(C9) 



(CIO) 



Substituting n q = n/q~ 9 / 2 we see, that indeed, the dom- 
inant part comes from the region of small qi. We can 
estimate all these integrals to get 



„ A 4 k 3 L 2 n 2 
2 

where we used the small q cutoff 1/L. Finally 



(Cll) 



S h k 3 L 2 n 2 I 



-'b 
Ik 



Pill 



(C12) 



and we conclude, that contribution from diagrams of type 
xG(ki + k 2 , uji + aj 2 )G(k + ki + k 2 , w + w\ + uj 2 ) (b) on Fig. 1(b) is much less, than contribution from one- 
xG(k + k, lj + uj 2 ) (CI) loop diagrams. But this is not the end of the story. Let 

us try to estimate contributions from diagrams of type 
where V a , Vb, V c , Vd are vertices, (e) on Fig. 1(b). Following the same guidelines, we obtain 



1 Remember, that in Zakharov-Sagdeev spectrum n(q) cx 
q~ 9 ^ 2 , and in Kadomtsev-Petviashvili spectrum n(q) tx g~ 4 . 
This assumption is true for the Zakharov-Sagdeev spectrum 
and is not true for Kadomtsev-Petviashvili one. 
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/71 1 
-^-A 4 khk 2 (k + ki)(k + h + k 2 )(k + k 2 )h kl h$\ 

[4] 
[5] 



X l(u k +U kl -UJk+k! +iTk,k 1 ,k+k 1 ) 
X (uJk + UJ kl + UJ k2 - LO k+kl+k2 + iTk l k 1 ,k 2 ,k+k 1 +k 2 

x (u>k + Lo k2 - u k+k , 2 + iTk,k 1 ,k+k 2 )] 1 



Let us again introduce k\ and k 2 as above, and substi- 
tuting h q = n/q~ 9 / 2 we obtain the following estimation: 



Ik 



7feC 2 



1 



(C13) 



Therefore we conclude, that contribution from two 
loop diagrams is dominated by planar diagrams and of 
the order of one loop diagrams contribution. 
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